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Abstract: An enabling debris capturing concept is experimentally demonstrated using 
an air-bearing testbed. To facilitate on-board implementation, the idea of artificial 
potential field is employed for the controller design of the debris capturing task. In 
particular, the artificial potential function is constructed such that all of its stable 
minima stratify the capturing condition, thereby avoiding the possibility of trapping 
into unwanted local minima. To experimentally simulate the mission of capturing 
space debris on ground, two air-bearing simulators are designed and developed, in 
which the first one is equipped with a single-degree-of-freedom robotic arm and the 
second one serves as a mock-up of the debris to be captured. The location and attitude 
of the simulators are measured by binocular cameras for control feedback. 


Experimental studies are performed to verify the results of the theoretic analysis. 
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1. Introduction 

Since the first space exploration mission in 1957, space activities have created an 
orbital debris environment, especially in the geo-stationary orbit, which severely 
threatens the safety of orbiting satellites and sustainable space development [1-4]. 
Some space activities have significantly produced the enormous amounts of orbital 
space debris which may remain in orbit for many years and interfere with future space 
activities[5]. Because of the nearly geometric growth in the amount of the space 
debris, many investigators are of the opinion that the debris density at the Low Earth 
Orbit (LEO) has already reached a sufficiently high level to trigger a cascading event 
(Kessler Syndrome) [6]. 

As a result, space debris has been of deep concern to the entire international 
space community during recent decades. In contrast to on-orbit servicing missions, 
almost all space debris is characterized by non-cooperative and tumbling motion, with 
no reflectors, markers or radio beacons, which would help with determining the 
relative position and attitude [7]. Therefore, three of the most challenging tasks in 
active debris removal (ADR) are minimizing the impact influence, de-tumbling and 
attitude synchronization, all of which have been studied widely in recent years [8]. To 
alleviate this situation, many capturing methods have been proposed. These methods 
can be classified into two main categories: contact and contactless capturing methods. 
Contact capturing, the most commonly proposed method, includes tentacles-based 


capturing [9], robotic arm capturing [10] and net capturing [11], etc. 


The robotic arm capturing method has attracted significant attention because of 
the advantages of stiff composite, higher Technology Readiness Level (TRL) 
technology and the capability of being easily tested on ground [8]. During the past 
decades, many kinds of the robotically capturing approaches have been proposed and 
investigated, however most of them mainly focus on theoretical and numerical studies 
[12]. Huang et al. established the impact dynamic model for target capturing and 
synthesized an adaptive control strategy to remove space debris by using a tethered 
space robot [13, 14]. Xu et al. proposed a method for autonomous rendezvous with a 
non-cooperative object in space, and developed a modeling and simulation system to 
verify the corresponding algorithms [15]. Lin et al. analyzed the polhode and 
tumbling frequency of the rotational space debris subjected to gravity-gradient torques 
and the orbit precession [16], which can be used to estimate the rotational state of the 
uniaxial space debris without any prior information except of the orbital parameters. 
Hu et al. proposed a trajectory planning algorithm capable of avoiding the 
configuration singularity of a space robot equipped with multiple manipulators [17]. 

Relatively less scholars have undertaken experimental research on this topic. For 
example, Dong and Zhu used a robotic arm installed on a fixed base to autonomously 
capture a mobile target with the help of visual servoing technique [18, 19], which 
required installing three non-collinear markers on the target at least. Moreover, many 
experimental studies using air-bearing spacecraft-simulator with robotic arm have 
also been reported. Toglia et al. demonstrated the control and sensing algorithms for 


assembly of spacecraft and space robots, based on the dual-arm two-link 


spacecraft-simulator system floating on a testbed [20]. Boning et al. built multiple 
satellite-manipulator system mounted on air-bearing microgravity simulators, to 
verify the construction of a large space structure by two or more space robots [21]. 
Virgili-Llop et al. demonstrated an onboard implementable optimization-based 
guidance approach for the capture of tumbling objects by spacecraft equipped with a 
robotic manipulator on air-bearing testbed [22]. 

A typical mission for space debris capture may consist of rendezvous, measuring, 
capturing and de-tumbling phases. In light of the restrictions due to laboratory 
environment, this study primarily focuses on the capturing phase. To simulate the 
mission of capturing space debris on ground, two air-bearing simulators are designed 
and developed, in which the first one is equipped with a single-degree-of-freedom 
robotic arm and the second one serves as a mock-up of the debris to be captured. A 
potential-energy based controller is proposed to complete an autonomous capture of 
the target mock-up that floats on an air-bearing testbed. The artificial potential 
function of the controller is constructed such that all of its stable minima stratify the 
capturing condition, thereby avoiding the possibility of trapping into unwanted local 
minima. The proposed control scheme was validated experimentally through a 
successful capture of the mock-up of the debris. 

The remainder of the paper is organized as follows. Section 2 describes the ADR 
problem of concern and the dynamic model of the ground-based experimental system. 
In Section 3, the design and stability analysis of the control law are presented. The 


experimental setup and results are given in Section 4. Finally, Section 5 draws the 


conclusions of the paper. 


2. Mission Concept and Dynamic Modeling 

Figure 1 give a schematic description of the conceptual mission for debris 
capture, in which a specially designed robotic arm is installed on a platform and 
employed to capture a space debris. The bottom surface of the end effector of the 
robotic arm is covered with numerous tiny bristles, like gecko’s foot, which could 
provide sufficient adhesion force for capturing debris. Previous studies showed that 
gecko foot hairs structure displayed excellent adhesion/release capacity even in the 
harsh space environment [23-25]. The bottom part (Disk A) of the end effector is 
connected via a revolute joint to the base part (Disk B) such that the relative rotation 
between the space debris and the end effector can be accommodated, where only 
single-axis rotation of space debris is considered for preliminary verification of 


conceptual feasibility. 
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Fig. 1. Illustration of capturing active debris using a robotic arm 
In this conceptual design, the space mission for active space debris capture 


includes the following four phases: the first step is the rendezvous phase in which the 


spacecraft approaches from a relatively faraway place to the target debris; the next is 
the measuring phase in which the rotation characteristics of the space debris are 
observed using the on-board measuring system; the third step is the capturing phase in 
which the target is captured by manipulating the robotic arm; the final step is the 
de-tumbling phase in which the relative rotation between the debris and the end 
effector is eliminated and then the whole system is de-tumbled or de-spun via thrust 
torque or other means. 

An air-bearing testbed system is developed in this work to experimentally 
demonstrate the basic concept of debris capture, as shown in Fig. 2. The experimental 
system is mainly composed of a binocular camera, an image processing computer, a 
granite platform and two air-bearing simulators (ABS) floating on the granite 
platform via air-bearings. The chaser simulator (CS) is equipped with a 
single-degree-of-freedom robotic arm and the target simulator (TS) serves as a 
mock-up of the debris to be captured. 
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Fig. 2. The ground-based experiment system 


As an initial proof of concept, the robotic arm in this work is simplified as a rod 
connected to CS through a single-degree-of-freedom revolute joint. The development 
of the adhesive end-effector is not considered in the represent work, and, for 
simplification, the free tip of the rod is equipped with a permanent magnet instead. A 
small piece of iron is fixed at the center of the top surface of TS, and the robotic arm 
of CS , as shown in Fig.3, is driven by a joint servo-motor and can be laid down to the 
horizontal plane for capturing TS at the position of the iron piece. With the 
consideration of the quite limited area of the granite table, TS is statically placed at a 
fixed position on the granite table and treated as a resident object to simplify the 


experimental implementation. 
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Fig. 3. The robotic arm driven by joint motor 
As shown in Fig. 4, two coordinate systems are employed for the dynamic 
modeling of the ground-based experimental system: the inertial frame OXY and the 


body-fixed frame 0,x,y, of CS. The attitude angle of CS with respect to the OX 


axis is denoted by 6. The radius of the ABSs is R and the length of the arm is 
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L+2R. The coordinate vector of point o, is denoted by Fr,= |X OY T for 
i=],2, where X, and Y, represent the X and Y-coordinates of the center of the top 
surface of the two ABSs, respectively. The angle between the robotic arm and the 
0,X% -axisis @ = 71/6, when the robotic arm is laid down to the horizontal plane. The 
solid line indicates the position vectors of the robotic arm when laid horizontally, and 
the black dots represent fluorescent markers for vision-based tracking. For the sake of 
clarity, introduce a virtual point C with the following absolute coordinates 
r, =[X,+(L+2R)cos(6,+@), Y,+(L+2R)sin(@, +a)" (1) 


which coincides with the position of the free tip of the robotic arm when laid 


horizontally. 


Fig. 4. Definitions of inertial and body-fixed reference frames 
To achieve the capture of TS, the following maneuvering phases were 
implemented, as shown in Fig. 5. In phases a and b, CS approaches to TS from the 
edge of the granite platform by adjusting its attitude and location in real-time. In 


phase c, CS arrives at the pre-capture state such that the position of the virtual point C 


coincides with that of the center of the top surface of TS. Finally, in phase d, the 


robotic arm is laid down to capture TS. 


a) Initial Approach Phase b) Adjustment Phase 


c) Pre-caputure Phase d) Capture Phase 


Fig. 5. The four phases of the capture mission on the air-bearing test bed 
Orbital motions and perturbations are neglected in the following analysis because 
the capture is completed in a relatively short period of time. The kinetic energy of CS 


can be written as follows 
T = sit + : JO (2) 
where m,and J, are the mass and the rotary inertia of CS, respectively. According 
to the Lagrange equation of the second kind, the dynamic equation of CS can be 
derived as 
Mq=u, (3) 


where M=diag|m, m J,| is the inertia matrix, q=[X, Y, 6] is the 


JT 
generalized coordinate vector, and u, =| F, F,, T,| is the generalized external 
force vector which consists of the forces acting on CS’s center in the X and Y 


direction and the moment acting on CS. 


3. Controller Design and Stability Analysis 

The control problem of current concern is to maneuver the chaser simulator such 
that the position of the virtual point C gets equal to that of the center of the top surface 
of TS. The pursued control objective is to guarantee, if the robotic arm is laid down to 
the horizontal plane, the free tip of the robotic arm can be exactly positioned at the 
center of the top surface of TS. 
3.1 The design of the controllers 

In the planar case, x=| X D0, Xo Ý, 6] denotes the system state of CS at 
time t. The desired trajectory of the center of the top surface of CS is represented by 
the dotted line circle in Fig. 4, which ensures that the distance between the centers of 
the top surface of CS and TS is equal to L+2R, the length of the robotic arm. The 
attitude angle 0, of CS is given by 

6, = 2arctan| Sy/(dx+ Jő +59") (4) 

where d=[6x, 6 yf is the relative position of two fluorescent markers on CS in the 
inertial frame OXY . Hence, 0 € (—n, T), i.e., there exists a singular point when 
O0 =-n or 0 =r. 

A potential-energy based control algorithm is devised based on the artificial 
potential field method to achieve the desired maneuvers of CS [26-28]. The artificial 


potential function, defined to adjust the position and attitude of CS, is defined as 
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follows 


V* = Bl(X,—X,)sin(@, +@)-(Y, -¥,)cos(6, + ay} +[r-(L+2R) | ©) 


where £ isa positive coefficient and r=./(X,-X,) +(¥,-Y,) . Accordingly, the 


potential-energy controllers based for position regulation can be formulated as 


P=, ae K,X, 
i ax, 
ov" 
Fy =-K, =——-K*ř 


where the coefficients K, and K, are positive coefficients. Besides, the 


potential-energy based controller for attitude regulation is designed as 
T =-K, — -Kð (7) 


where K, is also a positive coefficient. It is worth noting that the values 
ofsin(@,+@) and cos(6,+@) can be directly determined by using ôx and dy, 
thereby free of the singular problem of solving @ using Eq. (4). 
3.2 The stability analysis of the controllers 
To facilitate the analysis, one defines 
v =(X,—X,)sin(@, +a) —-(Y, -Y,) cos(@, + @) (8) 
v, =r—(L+2R) (9) 

Theorem 1. All of the stable minima of the potential function defined in Eq. (5) 
correspond to the capturing condition that v, =v, =0. 
Proof: The first-order derivatives of the potential function with respect to the three 


coordinates are set as zero to obtain the extremum, namely 


— =2fv,— +21, =0 (10) 


SP get nS = a1) 
E = (12) 
It is easy to verify that 
seek ~H) =52(X,-X,) (13) 
Fen ape) SK XS) (14) 


Thus, it is obvious that the solution of Eqs. (10)-(12) should satisfy either of the 


following two conditions: 


Condition A: 
v, =0 
15 
ee (15) 
Condition B: 
ov, ov, 
—Y,)-—(X,-xX,)=0 
OX, ( 1 5) ƏY, ( 1 2) T 
(L+2R)sin(6,+@) 
X-X, Pa ee — 


among which Condition B also implies r=(L+2R)/(f+1).Note that V* 20 all the 
time, thus, it is easy to conclude that V* arrives at its minimum if Eq. (15) is 
satisfied. It should be noted that Condition A defined in Eq. (15) implies that the 
position of the virtual point C coincides with that of the center of the top surface of 
TS. On the other hand, the Hessian matrix of V* is examined to determine whether 
the solution of Eq. (16) is the stable minimum of V*. The simplified Hessian matrix 


of V* is 
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(-28-1)cosg’+8+1 (-28-1)sin pcos g an 


H’=2| (-2B8-l1)singcosg (28+1)cosg’-£ -Bp sing (17) 
+ 


B+ 
Bp Bp Be 
rad are “(eH 


where 9=0,+7/6 and O=L+2R. The eigenvalues of the Hessian matrix 
are 28 ( p+h+2B +1) I(B +1) ,Oand 2£+2, respectively. Consequently, under 
the Condition B defined in Eq. (16), the solution of Eq. (13) is the unstable saddle 
point of V*. Thus, the stable minimum of V*ř is uniquely defined by Eq. (15), 
thereby avoiding the possibility of trapping into any unwanted local minima. 


Theorem 2. CS under the control laws (6) and (7) is asymptotically stable 


* 


Proof: Denote the solution of aa =0 as q,=[X,, Y, 0,Į , which satisfies the 
q 


dynamic model. A candidate of the Lyapunov function is constructed as follows 


L -276+ m (X? +Y?) ERY (18) 


y 


Obviously, L, 20 always holds. L,=0 if and only if V*=0 and the generalized 
velocity vector of CS is equal to the zero vector. Consider the first derivative with 


respect to time of L, 
L, =J,6,0,+m,X,X,+mYY,+ KV" (19) 
Combining the control laws (6) and (7), Eq. (19) can be simplified as 
L, =-K,X/-K,Y?-K,@ (20) 
which is always less than or equal to zero. Thus, the dynamics of CS is asymptotically 


stable under the control laws (6) and (7). 
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To evaluate the performance of the control scheme, numerical case studies are 
performed in MATLAB using the numerical integrator ode45. The length L was 
taken to be 0.3 m, and the mass, the rotational inertia and the radius of CS were set as 
9 kg, 0.05 kgm? and 0.1 m, respectively. The control parameters, K,, K,, and 
K, were set as 1, 4, 5, and 3, respectively. 

Note that the stable minima defined in Eq. (15) are not unique but all correspond 
to the capture condition and, consequently, the final location of CS depends on its 
initial condition. Hence, two different initial states of CS are used in case studies for 
comparison. The two initial states of CS were taken to be 
[ X,.%,,9,.X,,¥,,6, |=[0.2,0.3,0,0,0,0] and [0.3,0.3,0,0,0,0] , respectively. TS 
statically resides at the position [X,,Y,|]=[0,0]. Figure 6 depicts the simulation 
results for a period of 50 seconds. Clearly, the position and attitude of CS achieve a 
stable state in approximately 15 and 10 seconds, respectively. By checking the value 
of X,, Y, and @, itis found that the artificial potential function is practically equal 
to zero and the capture condition is satisfied, which implies that the free tip of the 
robotic arm is at the center of the top surface of TS when the robotic arm is laid down 


to the horizontal plane. 
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Fig. 6. Dynamic responses of CS subject to two different initial conditions 


4. Ground-based Experiments 
The proposed control algorithm is further validated experimentally based on the 
air-bearing testbed shown in Fig. 2. CS is equipped with a single-degree-of-freedom 
robotic arm and TS serves as a mock-up of the debris to be captured. Each satellite 
simulator is equipped with four normal and four tangential cold-gas thrusters to 
produce the normal forces and torque, respectively. Figure 7 shows the hardware 


layout of the ABSs. 
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Fig. 7. The hardware layout of the two air-bearing simulators 

The main part of the ABS weighs about 9 kg which includes about 1 kg of CO, 
in an on-board gas cylinder. The nominal pressure in the gas tank is about 
approximately 7 MPa. After flowing through the first stage decompression valve, CO2 
are divided into two-ways with the same pressure about 0.6 MPa. One pipe is directly 
connected to the three air bearings to allow the ABS to float above the granite 
platform, and the other pipe goes through the filter valve, which further reduces the 
pressure to about 0.2 MPa, to supply the eight nozzles, as shown in Fig. 8. Note that 
the electromagnetic valve used to control the nozzle switch is sensitive to the 
impurities in CO», therefore, the filter valve is necessary. In addition, the filter valve 
could ensure that the air supply pressures to the nozzles almost keep constant with 


different input pressures of the filter valve. For the given pressure 0.2 MPa, a normal 
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nozzle could provide a thruster force of about 0.05 N and a pair of tangential nozzles 
could generate the torque of about 0.008 Nm. 
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Nozzles Air Bearings 


Fig. 8. Schematic view of the on-board cold gas system 


As aforementioned, a single degree-of-freedom robotic arm is used in this paper. 
The robotic arm, which is driven by a servo motor, is a carbon fiber tube with a length 
of 50 cm, and its end effector is a permanent magnet. Because of the structural design, 
the robotic arm can only rotate in a fixed direction relative to CS. An on-board 
computer in the middle stage of the ABS and an on-board power system at the bottom 
of the ABS allow it to work without external dependencies. The ABS communicates 
with the image processing computer via a ZigBee module. Compared to the Wi-Fi 
module used by [29], the ZigBee module has lower power consumption and requires 
smaller battery capacity to achieve a longer usage time. As shown in Fig. 9, four 
fluorescent markers are placed on the top of the ABSs for the sake of optical 


measurement. 
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Fluoresent Markers 


Fig. 9. Fluorescent markers at the top of ABSs 

The positions of the fluorescent markers are captured in real-time by the 
binocular camera installed on the ceiling; subsequently, the acquired image data is 
further transmitted to the image processing computer. Whenever the computer obtains 
the three-dimensional coordinates by processing the digital image, it sends the 
information to the on-board computer mounted on the ABS via ZigBee wireless 
communication. According to the 3D coordinates of the markers, the feedback control 
law is solved by the on-board computer. By modulating the continuous control law by 
the Pulse Width/Pulse Frequency (PWPF) method [30], the switch command 
sequence for controlling the electromagnetic valves is obtained. Figure 10 shows the 
work flow chart of in the ground-based experimental system. It is worth noting that 
the single-degree-of-freedom robotic manipulator is also controlled by the on-board 
computer. Once the condition Ir. =r,| <0.02m is satisfied for more than 10 s, the 


robotic arm is put down to capture TS. 
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Fig. 10. The work flow chart of the experimental system 

The experimental test is performed over a period of 60 seconds. The 
experimental results shown in Figs. 11-14 are obtained using the binocular camera 
and the image processing system. TS resides statistically at the positions (0.238, 
—0.090) m in the inertial frame used for optical measurement. The control parameters 
K, K,, and K, were set as 24, 50, 2 and 1, respectively. 

Figures 11 and 12 depict the absolute coordinates of the virtual point C, where 
the dotted line represents the position of the center of TS. The coordinates of the 
virtual point C fluctuate slightly back and forth near the target value, due to the 
limited precision of the PWPF modulation of the thrust, as shown in the Figs. 11 and 
12. However, the range of the coordinate fluctuation of the virtual point C is so small 
and does not affect the capture task. It can be found from Fig. 13 that the value of v, 
converges to zero with slight fluctuations. As shown in Fig. 14, the relative distance 
between the centers of the chaser and target is close to 0.8 m at the initial moment, 
and gradually converges to the target value. After t=40s, the relative distance 


remains at approximately 0.5 m and the capture condition Ir. =T, | <0.02m is met. 
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Finally, the capture task can be completed by laying down the robotic arm to place its 


free tip at the ce 
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Fig. 11. The X coordinate of the virtual point C 
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Fig. 12. The Y coordinate of the virtual point C 
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Fig. 13. The time history of v, 
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Fig. 14. The relative distance between the centers of the two ASs 


5. Conclusions 


Numerical and experimental studies are performed to demonstrate an enabling 
debris capturing concept based on the robotic arm on-board a chaser spacecraft. The 


experimental system has a chaser simulator equipped with a robotic arm and a target 
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simulator as a mock-up of the debris. The robotic arm is simplified in this work as a 
rod connected to the chaser simulator through a single-degree-of-freedom revolute 
joint. The position and attitude controller of the chaser simulator is synthesized using 
the method of artificial potential field such that the free tip of the robotic arm can be 
exactly positioned at the center of the top surface of the debris simulator to achieve a 
successful capture. The artificial potential function is constructed to avoid the 
possibility of trapping into any unwanted local minima, and the Lyapunov method is 
applied to demonstrate the asymptotically stability of the control scheme. The results 
of the numerical and experimental studies well demonstrate the efficiency and 
performance of the proposed control scheme. The proposed concept and scheme are 
currently verified using a simplified model which only accounts for the 
two-dimensional system motions and a single degree-of-freedom robotic arm. 
Significant extensions are required to handle the challenges of more realistic scenarios, 
such as the relative motions in the three-dimensional space, the complex dynamics of 
multiple degree-of-freedom robotic arm with adhesive effector and the pose 


estimation of non-cooperative targets. 
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Highlights: 
1) The control scheme is synthesized using artificial potential field. 
2) The constructed artificial potential function is free of unwanted local minima. 


3) The proposed scheme is experimentally demonstrated using an air-bearing testbed. 


